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• ^ Individual contacts between people serve as pathways where infections may 

iP propagate. Characterizing and modeling these contacts is, therefore, funda- 

mental to unveil potential infection routes, to understand the emergence of 
epidemics, and to control or avoid the epidemics [BBVOS]. One way to repre- 
sent the patterns of contacts is through networks (BBV08[ INlOj . The network 
^ structure captures the heterogeneity of individual behavior and goes beyond 

random mixing where everyone can interact with everyone else. While static 
structures have been the common modeling paradigm, recent results suggest 
that temporal structures play different roles to regulate the spread of infec- 
tions [RLHlll ISKLlOi ISVBll) (or infection-like dynamics, as in information 
spreading [KKPllp . This has shifted the research attention to the under- 
standing of such changing structures (HS12] . On temporal networks a par- 
,—1 ticular vertex is active only at certain moments and inactive otherwise such 

that a contact is not continuously available [RLHIH IKKPIH IHS12| . The 
. ^ temporal order of contacts restricts the infection routes between vertices in 

comparison to non-temporal networks [TSMlOl IPSll] . The time between two 
^ consecutive vertex-activation events is not necessarily uniform, but in many 

datasets it follows heterogeneous patterns (e.g. bursts), in particular it does so 
on several networks of relevance to spread of infections, as for example, sex- 
ual |RLH10| , proximity [SKLIQI ISVBll] , or communication contact networks 
(electronic virus) [JE09, KKPll]. 

Previous research is mostly based on data-driven studies in which the evolu- 
tion of the epidemics is contrasted between the original temporal sequence and 
randomized versions, where the original time-stamps in the edges are shuffled 
to destroy the temporal correlations |RLH11[ IKKPll] . While this approach is 
adequate for many purposes, it is difhcult to study the sole effect of a particu- 
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lar temporal constraint. Other topological structures (e.g. degree distribution, 
degree-degree correlation or community structure) can be present and also con- 
tribute to regulate the infection spread. Some studies suggest that heteroge- 
neous inter-event time implies on a slow decay time of the prevalence during 
an infection dynamics jVRL07] and slowdown of the spreading dynamics in the 
context of communication networks |IE091lKKPll| . The opposite effect, i.e. the 
speed up of the infection growth due to broad inter-event activation times, is 
observed in a sample of sexual contacts network |RLHllj . To improve the under- 
standing of such contrasting behavior, in this chapter, we present a simple and 
intuitive stochastic model of a temporal network and investigate how epidemics 
co-evolves with the temporal structures, focusing on the growth dynamics of the 
epidemics. The model assumes no underlying topological structure and is only 
constrained by the time between two consecutive events of vertex activation, 
hereafter called vertex inter-event time. 

The network model consists on random activations of a vertex according to 
a pre-defined vertex inter-event time distribution. The vertices active at a given 
time are randomly connected in pairs during one time unit. The link is then 
destroyed and the vertices set to the inactive state. The infection event occurs 
through this link if one of the vertices is in an infective state. We study this 
model by using a susceptible infective dynamics with one (SI) and two (SII) 
infective stages. The first dynamics is motivated for being an upper limit case, 
where once infected, the vertex continues infecting at every contact |H00j . The 
second dynamics is more realistic and corresponds to a model of HIV spreading 
including an acute (high infectivity) and chronic (low infectivity) stages of in- 
fection with different periods [HAFOSj . If the second stage is set to zero in the 
SII model, we recover the susceptible infected recovered (SIR) dynamics [HOP] . 

To study our model, we compare the effects of two vertex inter-event time 
distributions on the spread of an infection. We consider the geometricj^ (corre- 
sponding to uniform probability of being active) and the power-law (to repro- 
duce heterogeneous inter-event patterns) distributions. The main observation is 
that the speed of the infection spread is different for both cases but the differ- 
ences depend on the stage of the epidemics. In comparison to the homogeneous 
scenario, the power law case results in a faster growth in the beginning but 
turns out to be slower after a certain time, taking several time steps to reach 
the whole network. 

The chapter is organized in 4 sections. We introduce the stochastic model 
and discuss some structural properties of the resulting network in Section [2] and 
introduce and describe the epidemics models in details in Section [3] The results 
of the SI epidemics on the networks following the different vertex inter-event 
time distributions are presented and discussed in Section |4] together with the 
results of the SII model. We conclude the chapter in Section [5] highlighting the 
major contributions. 



'^The geometric is a discrete time equivalent to the exponential distribution. 
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2 Network Model 



In this section, we describe the stochastic model of the temporal network and 
discuss some key structural properties of the evolving network. 

2.1 Network evolution 

The network is dynamic in the sense that vertices are active only at certain 
moments in time and inactive otherwise. As a consequence, edges also appear 
and disappear throughout the network evolution. In our model, whenever a 
vertex becomes active, it is randomly connected to another active vertex. The 
corresponding edge remains available during one time step, and is destroyed 
afterwards (see Algorithm IT]) . 



Algorithm 1 Network model 
Initialization step: 
for i = 1 ^ N do 

Tnexti ^ number drawn from -Dpow 
end for 

Network Evolution: 
for < = 1,2, ... do 

^next ^ ^next 1 

V ^ all nodes v such that Tnext^ = 

if IF I is odd then 

select one node u at random from V 

V ^ V\u 
end if 

while ^ do 

select two nodes at random u,v 
make a link between u and v 
Tnext^ ^ number drawn from -Dpow 
Tnext^ ^ number drawn from -Dpow 

V ^ V\{u,v} 
end vifhile 

end for 



The activation time is a fundamental property of our model and depends on 
the system of interest. It has been observed in different contexts that vertices 
are not active uniformly in time but often follows different waiting times, i.e. the 
inter-event time ATy between two vertex activations is not uniform but follows 
heterogeneous patterns. In particular, this inter-event time on sexual and prox- 
imity contacts, email and cell phone communication patterns are well described 
by power-law hke ATy distributions |ELHin[ ISKLTOl IS VBl II HEOfll IKKPllj . A 
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Figure 1: Probability of AT, the inter-event time between subsequent activa- 
tions, for the power-law with cutoff (full line) and the geometric cases (dashed 
line) on the a) vertices and b) edges. The axes are in log-scale. 

power-law ATy distribution means that there are bursts, i.e. trains of activations 
followed by periods of inactivity of various lengths. There are different theories 
trying to explain such behavior in the context of communication activity. One 
theory is the priority queuing model |B05) and the other is based on a nonho- 
mogeneous Poissonian process including periodic activity patterns |MSM08j . 

To simulate this characteristic in our model, for each active vertex at time 
t' , we sample the next activation time inoxt from a power-law (with cutoff) 
inter-event time distribution (eqn. [Tj hereafter referred to as -Dpow) such that 
ATy = tnoxt — t' ■ Note that this is equivalent of saying that the probability 
of being active at time tnoxt, given that the vertex was active at time i', de- 
creases as tnext gets larger and this decrease follows a power-law. Since we have 
a stochastic model, the cutoff is necessary to avoid that too large incxt's are se- 
lected during the evolution of the network. Since the epidemics dynamics occurs 
within few time steps, we set the cutoff to a reasonable value (Ai = 5 • 10^*) 
so that the exponential cutoff does not affect the dynamics and the power-law 
regime prevails. In line with empirical observations of the inter-event time (e.g. 
refs. |RLH10[ IKKPllj ). we set alpha to 2 (Fig. [T^). Random initial activation 
times <noxt are also selected from the distribution to all vertices, and to remove 
transient oscillations we let the network evolves during 1,000 time steps before 
the analysis (see Algorithm [T]) . 

P(AT) cx AT-"e"^i'^^, (1) 

We compare the effects of the heterogeneous inter-event time with a null 
model where the probability p of a vertex being active is constant (Poissonian 
process). This constant probability gives an exponential distribution of inter- 
event times (eqn. [2j hereafter referred to as Dg^o) where p ~ I/A2 (A2 is the 
mean inter-event time and set to be the same value as obtained from the I?pow 
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case, i.e. A2 — 4.024 for a = 2). To be more accurate in the simulations, we use 
a geometric distribution since we are working in discrete times. 

P(AT) cx e^^^^'^, (2) 

To generate random samples following the geometric distribution, one can 
use the inverse transform sampling, rounding down the value for the nearest 
integer |CSN09| . This function is typically available in math libraries for most 
programming languages. The case of the -Dpow distribution is a bit trickier and 
we use a method similar to the one proposed on ref. [ CSN09j , shortly described 
by Algorithm [2] 

Algorithm 2 Random values using a power-law with exponential cutoff 
vifhile X is rejected do 

r random number drawn uniformly from the interval [0, 1) 

X ^ LKin - ^) (1 - r)^i/("-i) + ij 
p e^^^, probability of acceptance of x 

p random number drawn uniformly from the interval [0, 1) 
ii p < p then 

return x 
else 

x is rejected 
end if 
end while 



We show on Figure [T^ numerically generated random samples using the 
methods described above for the two distributions I?pow and Dgco- It is visible 
that the exponential cutoff in the power-law has little or no effect for ATy values 
smaller than 300, but still the maximum values can reach at about ATy = 
10, 000, while the maximum values of the geometric distribution are not larger 
than 70. 

2.2 Characteristics of the temporal network 

The proposed model creates evolving networks without assuming any informa- 
tion about the network topology. After few time steps, the majority of the 
vertices have been in contact to every other vertex at least once, which means 
that the vertex-degree is approximatelly Ny. It is, therefore, more meaningful 
to highlight that the distribution of the total number of contacts has a charac- 
teristic value, although the spread is larger for -Dpow in comparison with the case 
of Dgeo- In other words, the vertices have a characteristic number of contacts 
with the same partner. The number of active vertices per time step nactivo is 
similar for both cases and for simulations using a total of 1, 000 vertices, ^active 
oscillates about a mean value of about 250 vertices (Fig. [2]) . 

The dynamics of the specific vertex inter-event times has different conse- 
quences on the edge inter-event time ATe, i-C- the time between the subsequent 
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Figure 2: Probability of having riactivo active vertices per time step in a network 
with a total of 1, 000 vertices. 

activation of the same edge. The power-law inter-event time on the vertices 
results in a distribution of ATe following a power-law functional form during 
a long interval, a characteristic also observed on empirical analysis of sexual 
contacts |RLH10| and cell phone communication |KKP11| . As expected, the 
uniform activation probability of the Dgco case gives a constant probability for 
the ATe distribution (Fig.[l)D). 

3 Epidemics models 

In this section, we describe epidemics models and present our simulation results 
for the spread of infections in the evolving network structure. 

3.1 SI and SII models 

To investigate the impact of the heterogeneous vertex inter-event time on the 
spread of infections, we simulate two complementary epidemics models on the 
evolving network. The first model, susceptible infective (SI), mimics a scenario 
where the vertices remain infective after an initial infection. Although this 
model is unrealistic to model actual infections, it is a good case study because 
it corresponds to an upper bound of infection or a worst-case scenario. To study 
the propagation of a hypothetical infection with a more realistic model, we run 
simulations for the susceptible infective model with two infective stages (SII). 
This model is adequate to represent the spread of HIV where an individual is 
highly infective during the initial months after a primary infection (acute stage), 
followed by a chronic stage of low infectivity. 

In both epidemics models, all vertices start the dynamics at a susceptible 
state, except for one random vertex in the infected state. A vertex becomes 
infective after a contact with another infected vertex with probability /3 (per- 
contact infection probability or transmissibility) . For the SI model, the states do 
not change after an event of infection, but for the SII, the probability of vertex 
A infecting vertex B is /3i in the initial Ti time steps after the first infection 
of A, and /32 after Ti time steps. Assuming that one time step corresponds to 



6 



one day, which is reasonable approximation for sexual contact networks, we set 
Ti = 90 days, which is the average time of the acute stage for HIV, and following 
typical values, we use /3i//32 = 10 jKJW971 iHAFOSj . We run the simulations 
on 100 network ensembles with 10 random infection seeds for each and thus we 
average our results over 1, 000 simulations. 



3.2 Mean-field SI epidemics for the geometric case 

In this section, we discuss the case of the temporal network where the vertex acti- 
vation follows a Poissonian process that can be described by a simple mean-field 
approximation. A discrete time SI epidemics evolving within a population of 
vertices homogeneously mixed can be described by the set of equations ([s]) |A94j . 



where St and It represent respectively the number of susceptible and infected 
vertices in the network at time i, N is the total number of vertices in the 
network (iV = St + It is constant for all times) , and 7 is the average number 
of infectious contacts an infected individual makes per unit time interval. The 
initial condition is given by /q > and Sq — N ~ Iq. 

The system of equation ([s]) describes the SI dynamics on our geometric case 
if we set AT — 1 and 7 — A/3, where A is the average number of contacts per 
vertex per time step, and /3 is the probability of per-contact infection between a 
susceptible and an infected vertex. Since each active vertex makes exactly one 
contact during a given time step, the average number of contacts per vertex per 
time step is, consequently, the probability for a given vertex to be active at a 
given time step (eqn. |4]) . 

A — Prob{a vertex is active at time t) = ^""^^'^"^ (4) 



4 Results 

In this section we present the results on the co-evolution of the epidemics model 
and the network structure, and discuss the impact of the heterogeneous vertex 
inter-event time on the spread of infections. 



4.1 Eff'ect of the inter-event time on the SI dynamics 

To obtain a global picture of the propagation of the infection, we show on 
Figure [3]i the average number of infected vertices per time, {il{t)). The figure 
contrasts the results for the different distributions of inter-event time Upow 
and I?geo- For all studied parameters, the infection grows faster in the initial 
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period for the case of Dpow in comparison to the case of Dgeo but slows down 
afterwards, such that it takes longer to infect all vertices in the case of a power- 
law distribution. The difference in the growth patterns is weaker for smaller 
values of /3. During the initial interval before the inflexion point, for both inter- 
event cases, the infection grows exponentially but at different rates. 




Figure 3: The average number of infected vertices, {n{t)), for 3 probabilities of 
infection, /3 = 1 (i), 0.1 (ii), and 0.01 (iii), for a) power-law Dpow and geometric 
Dgeo cases, and b) mean-field D^f and geometric -Dgco cases. The axes are in 
log-scale. 

Figure|3j3 shows the evolution of {^{t)) for the mean-field (Umf) approxima- 
tion and for the numerical simulation of the stochastic model using the geomet- 
ric vertex inter-event time. For the mean-field case, we set an initial condition 
Jo = 1 and riactivc — 249.75 (a value obtained from the numerical simulation, 
see Fig. pj) which gives A2 = ^ ^■'^^ ^ network with 1, 000 

vertices. J?he theoretical and numerical curves closely overlap for most of the 
interval of study, with a little mismatch when almost all the network is infected. 

The speed up in the early stage and the slowdown in the final stage of the 
infection growth for the power-law case in comparison to the geometric case is 
observed for other values of a, i.e. the exponent of the vertex inter-event time 
distribution (eqn. [l]) . Figure [4] shows the difference in the number of infected 
vertices in the power-law case, (r2(i)pow), and the number of infected vertices in 
the geometric case {il{t)geo) during the simulation of 1, 000 time steps. Positive 
values indicate that the power-law case infects more vertices than the geometric 
one, and negative values indicate the opposite behavior. Independently of the 
exponent, we observe that during a significant interval, the power-law case leads 
to a faster infection growth, and this interval increases with decreasing a. 

To better understand and characterize the differences in the growth in the 
early and final stages, we plot the number of time steps necessary to reach 
10% {TiQ%, Fig.[5]i) and 90% {Tgo%, Fig. (5)3) of the vertices. We compare how 
the mean and the median values of this measure vary with /3 for each vertex 
inter-event case. The median is a better statistics for non-symmetric broad dis- 
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Figure 4: Difference between the number of infected individuals if the vertex 
inter-event time fohows a power-law distribution and the number of infected 
individuals if the vertex inter-event time follows a geometric distribution, A5' = 
(f2(t)pow) — (^(i)gco)- The letter a corresponds to the exponent of the power- 
law, and the geometric version always uses the same average inter-event time 
as the power-law. The x-axis is in log-scale. 
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tributions, whereas both the median and the mean are adequate statistics for 
symmetric distributions. Since the distribution of Tio% has no characteristic 
value for the power-law case (see Fig. [5]: and discussion bellow) , we use both 
statistics to characterize the distributions and to compare the diverging effects. 
Figure 5^ shows that the mean and median for both D 
to similar values for small /3, and the median of Tio% for Dpow 
than for the case of -Dgco, indicating that the infection takes less time to reach 
10% of the vertices in the I?pow case. Conversely, due to the shape of the distri- 
bution (Fi g. [5p ) , the mean value gives a biased value and suggests the opposite 
effect (Fig. l5ti). To reach 90% of the vertices, however, the situation is different 
and the two statistics (the median and mean) are similar for both the case of 
-Dpow and for I?goo (Fig.jsja). At this later stage of infection, the behavior of the 
epidemics is more homogeneous and we identify a characteristic time in which 
90% of the vertices are infected. (Fig.jsji). These results indicate a slowdown on 
the infection growth in case of power-law inter-event time distribution (Fig.pp), 
a result that is in agreement with previous theoretical studies |VRL07l IIE09j . 

As discussed above, the I3pow case results in a broad distribution of times 
for early infections (Tio%), meaning that although there is a peak in Fig. [sj;, 
there are several possible scenarios of outbreaks, and these possible scenarios 
are mainly because of the heterogeneous vertex activity. In fact, out of 1, 000 
trials, 145 result in times above 200 time steps, and the largest one is 4, 453 
time steps. As a comparison, the largest time observed in the I?goo case is 42 
time steps. This narrow distribution of the Dgco cases shows that the choice of 
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Figure 5: The per-contact infection probability (3 versus the mean and median 
time for the infection to reach a) 10% and b) 90% of the network for both Dpow 
and Dgeo cases. In case of /3 = 1, the probability of c) 10% and d) 90% of the 
vertices being infected after, respectively, rio% and TgQ% time steps. The y-axes 
in a) and b) are in log-scale. 



the initial condition has little influence on the propagation of the infection. The 
large Tio% values observed in the Dpow case are responsible for the divergence 
between the mean and the median. A similar effect, i.e. the distribution of inter- 
event times for the I^pow case being broader than for the Dg^o case, is observed 
when 90% of the vertices are infected. This is not that strong however because 
for Tqq% the distribution is more symmetric in case of I?pow (Fig-[5|l)- Out of 
1,000 trials, 71 result in times above 800 time steps, being 4,843 the largest 
value. As a comparison, the largest value for the Z3gco case is 60 time steps. 

4.2 Effects of the inter-event time on the SII dynamics 

In the case of the SII dynamics, vertices go through two infective stages where 
the individual contribution to the global epidemics varies at each stage. Since 
we set the period of the first infective stage to Ti =90, the growth curve for 
Pi = 1 and =0.1 (Fig. [6| is not much affected in the early stage because 
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Figure 6: Average number of infected vertices {^{t)) in case of SII dynamics 
for different combinations of per-contact infection probabilities /3i and /32', (i) 
/3i = 1 and /32 = 0.1, (ii) /3i = 0.1 and /Sa = 0.01, and (iii) f3i = 0.01 and 
(32 = 0.001. The axes are in log-scale. 

most of the network is infected before the period Ti . More specifically, according 
to Fig. [5]d, the mean value to reach 90% of the network is at about 400 time 
steps. After the inflexion point, however, we see that the slowdown interval 
starts earlier in the SII epidemics. For smaller /3 the high and low infective 
regimes affect the infection growth significantly, especially, for f3i = 0.01 and 
/32 = 0.001. This effect is particularly interesting because it corresponds to the 
range of per-contact infection probability of HIV on several societies |KJW97j . 
For these values, the epidemics waits roughly 10 times longer to take off in 
comparison to a one infective stage scenario (T ^ 100 for SI, and T ^ 1, 000 for 
SII, see also Fig. [3^). 

5 Conclusion 

Recent research has suggested that temporal constraints on network structures 
affect the spread of infections. In particular, there is evidence that the time 
between two consecutive events of vertex activation, observed in some empirical 
networks, speeds up (or slows down in case of STlike models of information 
spreading) the spread of infections in comparison to homogeneous vertex inter- 
event times. These results are based on contrasting the infection growth on 
the original empirical network and on null models of networks where the origi- 
nal time-stamps are randomized by a reshuffling procedure. To investigate the 
changes in the shape of the infection growth more carefully, we propose in this 
chapter a simple and intuitive stochastic model of a temporal network where the 
topology is left completely random but the vertex inter-event time is explicitly 
controlled by using a power-law (with cutoff) and a geometric inter-event time 
distributions. 

By simulating a SI epidemics model co-evolving with the network structure, 
we observe that the inter-event time affects the infection growth differently 
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according to the stage of the epidemics. On early stages, in comparison to the 
homogeneous scenario (Poisson-hke dynamics), the heterogeneous (power-law) 
inter-event time results in a faster infection growth, while at later stages, the 
infection growth is slower. In our model, these differences results from the inter- 
event time and do not depend on the network topology (e.g. degree distribution, 
community structure, degree-degree correlations). We also observe that the 
heterogeneous inter-cvcnt time creates a diversity of outbreaks and highlights 
the influence of the initial conditions on the infection spread. In other words, 
while the homogeneous inter-event time results in characteristic times to reach a 
certain number of vertices, the heterogeneous produces several scenarios where 
the infection may take long times before affecting the network significantly. 
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